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摘 要 : 基于 若干 热 辐 射 输 运 求解 方法 的 一 维 辐射 平衡 和 温度 间断 算 例 结果 的 对 比分 析 , 提出 了 一 个 由 离散 坐标 法 、 
Pi 近似 和 光学 厚度 极限 近似 构建 的 热 辐 射 输 运 解 算 方 法 库 , 该 方法 库 可 与 再 入 飞行 器 热 化 学 非 平 衡 流 场 控制 方程 实 
现 有 效 耦 合 模拟 。 利 用 该 方法 库 模 拟 了 60km 高 空 、 马 赫 数 35 的 球 头 算 例 ， 计 算 表 明 热 辐射 对 高 温 流 场 具有 显著 的 
“冷却 ”作用 ; 辐射 传 热 与 对 流传 热 水 平 相 当 ， 流 场 非 耦 合 热 辐射 计算 的 辐射 热流 将 高 出 耦合 计算 值 27%， 此 类 飞 
行 条 件 下 的 数值 模拟 研究 应 当 考虑 热 辐射 耦合 效应 。 
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Abstract: Based on comparisons of the results of the one-dimensional radiation equilibrium and 
temperature discontinuity model cases, a thermal radiation solving method library consisting of the 
discrete ordinates method, P! approximation, and the optically limiting approximation is proposed. The 
library can achieve the effective coupled simulation with equations describing the thermo-chemical 
nonequilibrium flow over the reentry vehicle. In the 60 km altitude and Mach number 35 sphere model 
case, lt ls found that the thermal radiation can ‘coo!l’ the high temperature shock layer remarkably. It is 
also demonstrated that the radiative heating is comparable to the convective heating, and the radiative 
heating of the simulation uncoupled with the radiation overestimates 27% of the coupled value, which 
illustrates that coupling effects of the radiation should be considered in corresponding research. 
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] 一 维 切 平板 公式 求解 热 辐 


0 引 应 。Candler 和 ParkD 采 ) 
是 射 输 运 ， 灶 


再 入 飞行 器 是 载 人 登 月 等 深 空 探测 任务 的 主 


体 ， 其 以 极 高 速度 ( 约 10 km/s) 返 回 地 球 大 气 层 时 ， 


外 围 高 温 激 波 层 流 场 (可 达 10 K) 将 对 
酷 的 气动 加 热 ， 极 大 威胁 飞行 器 安全 


飞行 器 造成 严 
， 甚 至 导致 任 


务 失败 趾 。 因此, 准确 预测 再 入 飞行 器 气动 热 环 境 ， 


提供 热 防护 设计 依据 ， 已 成 为 相应 深 


空 探测 任务 所 


合 热 化 学 非 平衡 流 场 计生 
型 在 近地点 78 km 高 空 、9 knys 速 度 的 飞行 状态 下 驻 
点 辐射 传 热 约 6x10′ W/m?。Feldick 


给 出 AFE 模 


“分 别 利用 一 维 


切 平板 近似 、 蒙 特 卡 洛 模拟 和 Pi- 蒙 特 卡 洛 混合 方法 
与 DPLR(Data Parallel Line Relaxation) 程 序 结合 ， 对 
MPCV(Multi-Purpose Crew Vehicle) 返 回 舱 开展 非 平 


衡 流 场 耦合 热 辐射 模拟 ， 发 现 切 平 板 近似 预测 的 辐 


| 


临 的 关键 技术 难题 。 射 热流 比 蒙特 卡 洛 模拟 高 估 30%。 高 铁 锁 等 四 在 光 
高 速 再 入 飞行 器 高 温 大 气流 场 特性 极为 复杂 ， 学 厚度 薄 极 限 假设 下 沿 空 间 路 径 积 分 发 射 系数 获得 

分 子 振动 能 激发 和 多 组 分 化 学 反应 等 热 化 学 非 平衡 辐射 热流 ， 高 温 流 场 通过 求解 化 学 非 平 衡 N-S 方 程 
过 程 与 热 辐射 强烈 耦合 .此 时 飞行 器 气动 加 热 预 测 ， 得 到 ， 认 为 再 入 热 环境 严酷 区 中 返回 舱 辐 射 热流 贡 


除 考虑 传导 和 对 流 作用 外 ， 必 须 计 及 


热 重 射 耦合 效 


献 高 达 总 热流 的 30%。 


目前 ， 国 内 外 研究 多 采用 单一 的 已 有 热 辐 射 输 
运 求解 方法 直接 与 流 场 数 值 格式 开展 耦合 模拟 。 然 
而 ， 现 有 高 精度 方法 对 热 辐射 在 无 穷 维 光 谱 和 传播 
方向 上 的 刻画 ， 无 疑 会 给 耦合 模拟 带 来 庞大 的 运算 
量 和 存储 量 ， 而 简化 热 辐 射 模型 又 存在 预测 精度 不 
足 的 问题 。 因 此 ， 单 一 热 辐射 输 运 求解 方法 无 法 满 
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中 当前 再 入 飞行 器 气动 热 环境 高 效 模拟 需求 。 本 文 
基于 一 维 热 辐射 输 运 模型 算 例 分 析 ， 充 分 考虑 再 入 
飞行 条 件 下 流 场 的 光学 厚度 和 数值 格式 特征 ， 构 建 
了 适用 与 再 入 飞行 器 热 化 学 非 平衡 流 场 方程 耦合 模 
拟 的 热 辐射 输 运 解 算 方法 库 ， 并 利用 该 方法 库 成 功 
实现 了 Ma =35 再 入 飞行 条 件 下 球 头 算 例 的 流 场 热 
辐射 耦合 模拟 。 


1 一 维 热 辐射 输 运 模型 


采用 统一 的 一 维 模 型 算 例 对 现 有 若干 热 辐射 输 
运 求解 方法 开展 对 比分 析 ， 模 型 算 例 几何 描述 同文 
献 [5]， 示 意见 图 1: 两 无 穷 大 各 向 同性 漫 射 平行 壁 
面 间 充 满 参与 性 灰 介 质 ， 介 质 热力 学 和 辐射 特性 与 
温度 无 关 。 平 行 壁面 间距 二 ， 法 向 坐标 ， 光 学 厚度 
坐标 z =xy，x 为 介质 辐射 吸收 系数 。 壁 面 1、2 的 温 
度 分 别 为 Tu 、7Tia， 发 射 率 分别 为 sf、s 。 辐 射 强度 
/与 y 轴 方向 夹 角 y， 周 向 角 pg。 该 模型 具有 一 维 轴 对 
称 性 质 ， 不 妨 取 等 间隔 空间 离散 ， 离 散 点 总 数目 
N+1， 物 理 间 隔 Ay =L/N， 光 学 厚度 间隔 At =tj/N。 


2 各 方法 的 一 维 模型 算 例 对 比分 析 


选取 目前 研究 中 常用 的 多 种 热 辐射 输 运 求解 方 
法 进行 一 维 模 型 算 例 对 比分 析 研 究 。 本 文 考察 的 方 
法 有 切 平板 近似 (Tangent Slab Approximatiom)、 蒙 特 
卡 洛 模拟 (Monte Carlo Method) 、 离 散 坐 标 法 
(Discrete Ordinates Method) 、 Pi 近似 (Pi 
Approximation) 以 及 两 种 光学 厚度 极限 近似 方法 : 光 
学 薄 近 似 (Optically Thin Approximation) 和 光学 厚 近 
似 (Optically Thick Approximation)。 切 平板 近似 实际 
上 是 热 辐射 输 运 方程 的 一 维 积分 解析 解 ， 对 于 图 1 
所 示 情 形 ， 切 平板 近似 有 辐射 热流 : 


gq, (7)=2[ TE,(7)-JE,(r, -7)| 


42c| fT (OB, (rn) a TO) E(t-7) | () 


其 中 7 表示 温度 , ;和 Es 分 别 为 二 阶 和 三 阶 且 指数 积 
分 函数 外， 且 有 


蒙特 卡 洛 模拟 结果 在 样本 量 足 够 大 时 亦 可 认为 是 热 
辐射 输 运 的 准确 解 加， 因此 本 文 一 维 算 例 模拟 可 以 
切 平板 近似 和 蒙特 卡 洛 模拟 的 结果 作 基 准 。 光 学 薄 
近似 一 般 要 求 特征 光学 厚度 tz <<1， 而 光学 厚 近 似 
则 要 求 z >>1。 各 方法 的 具体 原理 和 详细 表达 可 参 
见 文献 [6] 和 [7]。 本 文通 过 FORTRAN 自 编程 实现 上 
述 各 方法 计算 。 
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图 1 一 维 热 辐射 输 运 模型 几何 示意 图 
Fig.1 Sketch of one-dimensional radiative transfer model 

2.1 辐射 平衡 

假定 图 1 中 两 壁面 间 介 质 辐 射 平 衡 , 辐射 热流 9 
空间 梯度 为 零 ， 即 dg/dy =0 W/m ; 取 壁 面 温度 Ti 
=1000 KK，7Tw =700 K， 发 射 率 e1 =e, = 1.0， 分 别 计 
算 x =0.7、1.5、3.0 mo 的 工 涡 。 取 L =1.0 m。 定 义 无 
量 纲 辐射 热流 : 


A 


a g; 
本 _ 


为 确定 离散 点 数目 N， 针 对 基准 方法 切 平板 近似 开 
展 网 格 无 关 性 计算 ， 结 果 如 图 2 所 示 。 图 2 说 明 ， 
离散 网 格 点 数目 >30 时 ， 切 平板 近似 的 辐射 热流 
值 基本 稳定 ， 同 时 无 量 纲 温度 分 布 数据 也 不 再 随 网 
格 加 密 而 变化 。 因 此 ， 本 文 一 维 算 例 的 取 值 均 大 
于 30， 以 保证 切 平板 近似 方法 计算 的 准确 性 。 本 节 
辐射 平衡 算 例 取 N=50。 
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Fig.2 Grid independent simulation for tangent slab 
approximation: (a) nondimensional radiative heat transfer, (b) 
distribution of nondimensional temperature 

图 3 展示 了 本 节 辐 财 平 衡 算 例 各 方法 的 温度 分 
布 结果 。 离 散 坐 标 法 (DOM) 与 基准 方法 切 平板 近 似 
(TSA) 以 及 蒙特 卡 洛 模 拟 (MC) 的 温度 值 几 乎 一 致 ， 
Pi; 近似 (Pi) 计 算 值 也 与 基准 解 接近 。 光 学 薄 近 似 
(THIN) 的 温度 为 固定 值 ， 而 光学 厚 近 似 (THICK) 同 
切 平板 近似 的 数据 也 存在 很 大 差异 ， 但 随 着 总 光学 
厚度 值 专 增 大 ， 光 学 厚 近似 预测 效果 会 逐渐 改善 。 
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3 一 维 辐射 平衡 问题 各 方法 温度 分 布 结果 


Fig.3 Temperature distributions of one-dimensional radiative 


equilibrium calculated by different methods 


各 方法 的 无 量 纲 辐射 热流 


差 为 以 切 平板 近似 热流 值 为 基准 计算 的 相对 误差。 
切 平板 近似 与 蒙特 卡 洛 模拟 两 类 基准 方法 的 差异 极 


小 ， 离 散 坐 标 法 对 辐射 热流 预测 类 


甫 度 很 高 ， 而 Pi 


近似 的 辐射 热流 误差 均 在 5% 以内， 但 两 种 光学 厚 
度 极 限 近 似 的 热流 值 误 差 均 很 大 。 


表 1 
Table 1 


Method t=0.7 


TSA 0.634 
MC 0.634 
DOM 0.635 
Pi 0.656 
THIN 1.000 
THICK 1.905 


各 方法 无 量 纲 平衡 辐射 热流 值 对 比 


Comparison of non-dimensional equilibrium 


radiative heat fluxes calculated by different methods 


Error t=1.5 Error t=3.0 Error 


0.000E+00 0.457 0.000E+00 0.301 0.000E+00 


8.883E-06 0.457 1.734E-04 0.302 6.968E-04 


7.150E-04 0.458 1.304E-03 0.302 2.978E-03 


3.369E-02 0.471 2.915SE-02 0.308 2.077E-02 


5.764E-01 1.000 1.187E+00 1.000 2.318E+00 


2.003E+00 0.889 9.440E-01 0.444 4.744E-01 


2.2 温度 间断 


仿照 再 入 飞行 器 高 温 激 波 层 流 场 ， 设 计 温度 间 


断 算 例 : 


两 避 间 介 质 温 度 分 为 均匀 两 层 ， 分 段位 


置 y/L =0.3， 


同 为 10000 K; 近 壁 


近 壁 面 1 处 为 高 温 层 ， 温度 与 辟 温 也， 
2 处 为 低温 层 ， 温 度 与 壁 温 


Tw 均 为 300 KK。 分 别 计算 x=0.01、 0.1、 1.0、5.0m 


的 工 况 。 取 工 =1.0m，N=1000。 定 义 无 量 纲 辐射 能 
量 源 项 : 

dgq,. /dy 

7 4 

xo (74 -7%) Wy 

各 方法 无 量 纲 辐射 能 量 源 项 计算 值 见 图 4。 所 

有 光学 厚度 条 件 下 ， 离 散 坐标 法 与 基准 方法 切 平 板 

近似 结果 高 度 一 致 ，P; 近似 稍 有 偏差 。 小 光学 厚度 


t=0.01 时 ， 光 学 薄 近 似 与 切 平板 近似 结果 相近 ， 
但 随 去 增 大 ， 逐渐 失去 合理 性 。 光 学 厚 近 似 辐射 能 
量 源 项 值 始终 为 零 。 
3.0 c 
~ 2.0 E 
3 [ 
号 10 上 
.2 [ 
中 00G8 E6060 
B10 EF 一 重 - 一 一 TSA 
E 四 一 一 全 一 一 DOM 
[=| 上 一 一 会 一 Pi 
-2.0 —@— THIN 
EF ——E-—. THICK 
-230 人 
0.000 0.002 0.004 0.006 0.008 0.010 
tT 
(a) Tt =0.01 
3.0 5 
、2.0 上 
3 9 [ 
号 10 上 
.2 
中 008 E6000600GHE06000 
£ EF 
Ef 和 
只 -2.0 ee — i 
EF ——O—- THICK 
0 人 


MD 


nondimensional mw 
Dea] 


0.00 0.02 0.04 0.06 0.08 0.10 
(b) zt =0.1 

0 

中 

.0 上 

0 和 eee-e EEOC 


nondimensional oO 


-了 P， 


TSA 
DOM 


THIN 
THICK 


-3.0 - 
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 
T 


(d) t=5.0 


图 4 温度 间断 各 方法 无 量 


量 纲 辐射 能 量 源 


项 结果 


Fig.4 Non-dimensional radiative energy source terms of 


temperature discontinuity using different methods 


ti =0.1 时 各 方法 的 无 量 


表 中 误差 仍 是 以 切 ; 
对 误差 。 离 散 坐标 
相同 ，Pi 近似 也 能 给 出 正确 
件 ， 光 学 薄 近 似 结果 具有 一 定 可 靠 性 ， 但 光学 厚 近 
b 零 热流 结果 。 其 他 总 光学 厚度 条 件 
结果 表现 出 类 似 特 征 。 

tL =0. 1 各 方法 两 壁面 处 无 量 纲 辐射 热 


似 只 能 给 
计 复 
表 2 


纲 辐射 热 
“ 板 近似 热流 值 为 3 
法 与 切 平板 近似 的 辐 
值 。 此 类 小 光学 


% 沅 


基准 计生 
晶 射 热流 基本 


值 见 表 2， 
的 相 


厚度 条 


下 ， 


流 结果 对 比 


Table 2 Non-dimensional radiative heat transfers at 


two walls calculated by different methods with z, =0. 1 


Method 
TSA 
DOM 
Pi 
THIN 
THICK 


wall 1 Error 

0.944 C0.00E+00 
0.943 7.38E-04 
0.963 2.05E-02 
1.000 5.93E-02 
0.000 1.00E+00 


wall 2 
0.878 
0.876 
0.892 
1.000 
0.000 


3 热 辐 射 输 运 解 算 方法 库 


合 模拟 需求 站。 


方法 构建 一 


一 个 


本 文 提 


热 辐 射 输 运 解 算 方法 库 ， 


:能 差异 显著 ， ee 
出 如 下 解决 思路 : 集成 多 入 


需求 并 分 阶段 地 选择 库 中 方法 执行 热 帮 


Error 
0.00E+00 
1.45E-03 
1.58E-02 
1.39E-01 
1.00E+00 


Ne edt rh ei 合 热 辐射 求解 
必然 耗费 更 大 的 计算 资 
输 运 求解 方法 忆 


[各 类 热 辐射 


依 特 性 、 按 


下 射 耦合 模 


拟 。 解 算 库 4 


Si 


可 诊 


方法 应 满足 : 


E 确 预测 辐射 热 


成 本 ， 能 与 流 场 数值 格式 协同 求解 ， 


合 大 规模 并 行 计算 ,基于 2 节 模 型 算 例 对 比分 析 ， 


虑 各 方 法 程 ) 万 
性 能 特点 列 于 表 3， 表 : 


量 源 项 与 辐 
切 


S| 


虽 射 热 


F 板 近似 是 一 维 


求解 执行 过 程 ， 


总 结 出 


EN 二 口 


\ 流 计算 而 言 。 


预测 精度 主要 针对 辐射 能 


让 分 解析 解 ， 针 对 一 维 问题 


求解 是 精确 的 ， 但 其 只 能 沿 一 个 维度 路 径 积分 ， 攻 


而 用 于 高 维 计算 有 失 合 理性 ， 精 度 难 以 界定 。 蒙 特 


卡 洛 模拟 在 样本 量 充足 时 可 认为 精确 解 ， 但 会 相应 


带 来 计算 量 的 问题 。 离 散 坐标 法 在 一 维 求解 时 与 基 
准 方法 切 平 板 近似 和 蒙特 卡 洛 模拟 结果 高 度 一 致 ， 
是 基于 完整 的 辐射 传输 方程 的 直接 离散 方法 ， 且 容 


易 实 现 高 维 扩展 。P; 近似 求解 的 是 简化 的 辐射 传输 
方程 ， 精 度 适中 。 光 学 薄 和 光学 厚 近似 只 有 在 极 小 


(一 般 z<0.01) 或 极 大 光学 厚度 条 件 下 才 合 理 。 


表 3 


各 类 热 辐 射 输 运 求解 方法 特点 对 比 


Table 3 Characteristics of different methods solving 


radiative transfer 


方法 预测 精度 存储 量 计算 量 级 高 维 扩展 ”并 行 化 
可 
TSA 仅 一 维 精确 小 NY 复杂 
需 大 量 传 值 
MCM 精确 很 大 “一般 很 大 ”容易 可 行 
容易 
DOM 高 适中 NxM 容 
传 值 量 适中 
容易 
Pi 适中 小 N 容 
传 值 量 小 
THIN ”接近 光学 厚度 很 小 容易 
小 容易 
THICK 极限 条 件 时 合理 无 需 循环 传 值 量 很 小 


在 计算 资源 占用 方面 ， 切 平板 近似 沿 一 维 路 径 
进行 数值 积分 ， 存 储量 不 大 ， 但 在 各 网 格 点 处 均 要 
进行 上 、 下 游 逐 点 积分 , 一 轮 计算 需 完成 Y 次 代数 


运算 ， 其 中 为 离散 网 格 点 数 ， 计 算 量 相当 可 观 。 
蒙特 卡 洛 模拟 需要 追踪 每 个 能 束 吸 收 和 发 射 的 历 


程 ， 往 往 巨 大 的 能 束 样本 量 会 占用 很 大 的 存储 资源 


和 计算 量 。 离 散 坐标 法 需要 存储 各 角 方 向 上 的 辐射 


强度 , 在 每 个 网 格 点 处 需求 解 离散 角 方 向 数目 M 个 
代数 方程 ， 需 占用 一 定 的 存储 和 计算 成 本 。P; 近似 
仅 求解 关于 投射 辐射 (incident radiatiom) 的 二 阶 微分 


方程 ， 存 储 和 计算 量 均 不 大 。 光 学 厚度 极限 近似 由 
显 式 计算 辐射 量 ， 求 解 成 本 很 低 。 表 4 


温度 场 直 接 


列 出 了 求解 


2 节 中 辐射 平衡 算 例 各 方法 所 消耗 的 


CPU 时 间 ， 


其 中 蒙特 卡 洛 模拟 所 需 的 能 束 样 本 量 为 


100 万 个 ， 离 散 坐标 法 考虑 的 角 方 向 为 10 个 。 从 表 
4 中 可 以 看 出 ,各 方法 所 耗费 的 实际 计算 时 间 与 表 3 


的 定性 分 析 是 一 致 的 。 切 平板 近似 和 蒙特 卡 洛 会 


消耗 大 量 计算 时 间 , 而 离散 坐标 法 和 Pi 近似 计算 效 


率 较 高 ， 而 光学 厚度 极限 方法 几乎 不 消耗 时 间 。 
除 切 平板 近似 难以 推广 至 高 维 外 ， 其 余 方法 则 


易于 扩展 至 高 维 情形 。 同 时 ， 切 平板 近似 还 会 在 并 


行 求解 过 程 时 出 现 跨 区 域 积分 问题 ， 需 要 各 3 


行 区 


域 间 进 行 大 量 传 值 。 离散 坐标 法 和 Pj 近似 均 是 传递 
相 邻 区 域 交 界 处 的 热 辐射 相关 量 。 光 学 厚度 极限 近 
似 则 只 需 传 递 相 邻 区 域 交界 处 温度 值 。 

表 4 辐射 平衡 算 例 各 方法 所 消耗 CPU 时 间 (单位 : s) 


Table 4 CPU time consumption of solving radiative 


equilibrium case by different methods (unit: s) 


综 上 ， 


Method CPU time/s 
TSA 8.7813 
MC 10.7813 
DOM 0.3906 
Pp 0.1719 
THIN 0.0000 
THICK 0.0000 


本 文 放弃 不 利于 高 维 求解 和 并 行 处 理 的 
切 平板 近似 与 高 计算 成 本 的 蒙特 11 


FE 洛 模拟 。 本 文 热 


辐射 输 运 解 算 方 法 库 最 终 选 定 由 离散 坐标 法 、P1 近 
似 和 光学 厚度 极限 近似 ( 含 光学 薄 和 光学 厚 近 似 ) 三 


类 不 同 层次 的 方法 构成 。 该 方法 库 面 对 热 辐 射 输 运 
求解 的 一 般 执行 原则 为 ， 与 流 场 玫 
解 时 直接 选用 离散 坐标 法 ; 与 流 场 进行 耦合 求解 时 ， 
6 学 厚度 满足 极限 条 件 ， 可 先 用 光学 厚度 
合 模拟 ， 收 敛 后 转 用 离散 坐标 法 以 
节约 计算 量 ， 若 不 满足 光学 厚度 极限 条 件 ， 则 先 用 


若 流 场 的 》 


极限 近似 进行 夺 


F 展 非 耦合 


六 


者 确 求 


Pi1 近似 耦合 求解 ， 收 敛 后 转 用 离散 坐标 法 。 上 述 使 
用 原则 如 图 $ 所 示 。 
然 己 遇 答 运 解 代 方 法 床 
入 的 计算 量 
精 
确 
计 
算 
流 场 计算 收 人 
直 
接 
光学 砷 应 极限 近似 总 近似 
YY 
eT 痪 改 生生 尖 
计生 收 人 
给 出 辐射 结果 
图 5 热 辐射 输 运 解 算 方 法 库 操作 原则 


Fig. 3 Operation procedures of the thermal radiation solving 


method library 


4 再 入 飞行 算 例 粳 合 模拟 演示 


选取 Ma =35 的 R=1.0 m 球 头 模型 开展 再 入 飞 
行 条 件 下 高 温 流 场 与 热 辐 射 耦合 模拟 ， 演 示 本 文 构 


造 的 热 辐射 输 运 解 算 方 法 库 与 流 场 数值 格式 耦合 求 
解 的 有 效 性 。 取 自由 来 流 密度 pw =3.1x10“kgm， 
温度 思 =247 K， 大 气质 量 分 数组 成 p=0.77，Yo， 
=0.23， 上 述 取 值 对 应 60 km 高 空地 球 大 气 参 数 。 高 
温 流 场 求解 采用 文献 [8-10] 中 给 出 的 数值 格式 ,采用 
Park 两 温度 热力 学 非 平衡 模型 5 和 Gupta 的 11 组 
分 20 反应 化 学 反应 动力 学 模型 中 。 热 辐射 特性 计 
算 采 用 Anderson 两 台阶 模型 中 , 壁 温 取 ,=1000 K,， 
壁面 发 射 率 。 =1.0， 取 完全 催化 壁面 条 件 。 流 场 与 
热 辐射 而 合算 法 见 文献 [8]。 流 场 求解 和 热 辐 射 解 算 
方法 库 均 采 用 FORTRAN 自 编程 实现 。 计 算 网 格 量 
60x70， 按 照 气动 热 预 测 需求 ， 壁 面 法 向 第 一 层 网 
格 间距 1.0x10 mb3l。 该 球 头 模型 计算 量 相对 较 小 ， 
因此 流 场 耦 合 热 辐 射 求解 可 直接 采用 解 算 方 法 库 中 
高 精度 的 离散 坐标 法 开展 。 

图 6 给 出 了 流 场 耦合 热 辐射 计算 的 压强 和 温度 
分 布 ， 此 时 激 波 层 很 薄 ， 激 波 脱 体 距离 不 到 球 头 半 
径 的 1/15, 波 后 驻 点 区 域 平 动 - 转 动 温度 和 振动 ~- 电 
子 温度 均 超过 10000 K。 图 7 展示 了 流 场 内 热 辐射 
量 分 布 结果 ， 辐 射 能 量 源 项 和 平均 吸收 系数 分 布 形 
式 完全 对 应 于 流 场 的 温度 分 布 。 投 射 辐射 最 大 值 虽 
然 分 布 在 激 波 层 内 ， 但 高 投射 辐射 值 分 布 也 扩展 到 
了 激 波 上 游 来 流 ， 更 具 扩 散 性 。 

Ma =35 飞行 状态 下 , 辐射 能 量 源 项 已 具有 相当 
量 级 ( 约 105 WAmm2)， 热 辐射 对 流 场 作 用 显著 ， 但 激 
波 层 的 光学 厚度 仍然 很 小 , 不 足 0.005。 因 此 ， 本 文 
还 采用 了 解 算 方法 库 中 的 光学 薄 近 似 与 流 场 开 展厅 
合计 算 。 图 8 比较 了 流 场 非 看 合 热 辐射 (uncoupled)、 
耦合 光学 薄 近 似 (coupled thim 和 耦合 离散 坐标 法 
(coupled_dom) 的 驻 点 线 温 度 结果 。 与 非 耦 合 相 比 ， 
耦合 计算 的 激 波 层 温 度 更 低 ， 激 波 层 更 薄 ， 凸 显 了 
辐射 冷却 作用 四。 光学 薄 近 似 比 离散 坐标 法 而 合计 
算 的 辐射 冷却 效果 更 明显 ， 原 因 在 于 光学 薄 近 似 忽 
咯 了 辐射 吸收 作用 ， 仪 考虑 辐射 发 射 ， 辐 射 能 量 源 
项 值 更 大 ， 对 流 场 作用 更 显著 。 


中 | 


了 
iV 会 
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近似 和 离散 坐标 法 计算 的 驻 点 线 热 辐射 量 。 耦 合计 
算 的 热 重 射 量 绝对 值 均 小 于 非 耦合 值 ， 光 学 薄 极 限 
近似 耦合 计算 的 热 辐射 量 绝对 值 减 弱 更 明显 。 
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(b) (c) 
图 6 ”耦合 热 辐 射 的 流 场 模拟 结果 : (a) 压强 (单位 : Pa), (b) 
平 动 -转动 温度 (单位 : K), (c) 振动 -电子 温度 (单位 : K) 


Fig.6 Flowfield simulation results coupled with radiation: (a) 


一 


pressure (unit: Pa), (b) translational-rotational temperature (unit: 
K), (c) vibrational-electronic temperature (unit: K) 
表 5 给 出 了 流 场 非 耘 合 与 耦合 光学 薄 近 似 和 离 
散 坐 标 法 计算 给 出 的 驻 点 对 流 热流 (qe) 和 辐射 热流 
(g)) 值 。 非 看 合 辐射 热流 值 是 在 流 场 的 收敛 结果 基础 
上 ， 再 以 离散 坐标 法 进行 求解 得 到 的 ， 光 学 注 近 似 
无 法 给 出 辐射 热流 值 ， 其 辐射 热流 同样 是 在 耦合 计 


图 9 分 别 给 出 了 流 场 非 耦合 热 鲍 射 与 耦合 光学 
薄 近 似 和 耦 合 离散 坐标 法 计算 的 驻 点 线 N 和 e 质 
量 分 数 分 布 。 耦 合计 算 给 出 的 N 和 的 质量 分 数 明 


算 收 敛 后 的 流 场 基础 上 ， 以 离散 坐标 法 求解 获得 。 
由 表 5 可 知 ， 非 耦合 与 耦合 热 重 射 给 出 的 对 流 热 流 
直接 近 , 但 非 耦 合 给 出 的 辐射 热流 明显 高 于 耦合 值 ， 


显 低 于 非 耦 合 值 ， 原 因 可 能 是 辐射 冷却 作用 抑制 
电离 反应 。 因 此 ， 高 速 再 入 飞行 条 件 下 高 温 大 气 化 
学 反应 和 热 辐射 存在 显著 的 耦合 效应 。 还 可 以 看 到 ， 
光学 薄 近 似 耦 合计 算 辐 射 冷 却 效果 更 显著 ， 其 N 
和 e 质量 分 数 低 于 离散 坐标 法 耦合 计算 值 。 

图 10 给 出 了 流 场 非 耦合 热 辐射 与 耦合 光学 薄 


高 出 光学 薄 近 似 耦 合 结果 48%， 高 出 离散 坐标 法 耦 
合 值 27%， 光 学 薄 近 似 耦 合计 算 较 离散 坐标 法 耦合 
计算 的 辐射 热流 值 低 。 此 时 ， 流 场 的 辐射 传 热 可 与 
对 流传 热 水 平 相当 。 因 此 ， 开 展 此 类 再 入 飞 条 件 下 
气动 热 研究 时 ， 必 须 考 虑 流 场 的 热 重 射 耦合 效应 。 
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7 流 场 中 热 辐 射 量 分 布 : (a) 辐射 能 量 源 项 (单位 : W/m’)， 
(b) 平均 吸收 系数 (单位 : m 0, (c) 投射 辐射 (单位 ，W/m?) 


Fig.7 Contours ofradiation variables of coupled radiation 


simulation: (a) radiation source (unit: W/m’), (b) total 


absorption coefficient (unit: m !), (c) incident radiation (unit: 


W/m’) 
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8 ” 流 场 非 而 合 与 而 合 热 辐射 计算 的 驻 点 线 温度 


Fig.8 Temperature distributions along the stagnation line of 
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uncoupled and coupled radiation simulations 
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Fig.9 Mass fraction distributions along the stagnation line of 


uncoupled and coupled radiation simulations: (a) N', (b)e” 
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投射 辐射 (单位 ，W/m”) 


Fig. 10 Radiation results along the stagnation line of 


uncoupled and coupled radiation simulations: (a) radiation 
source (unit: W/m’), (b) total absorption coefficient (unit: m '), 


(c) incident radiation (unit: W/m’) 


9 流 场 非 耦合 与 耦合 热 辐射 计算 驻 点 线 组 分 质量 分 数 结 


10” 流 场 非 契 合 与 看 合 热 辐射 计算 驻 点 线 热 辐射 结果 : (a) 
辐射 能 量 源 项 (单位 : W/m’), (b) 平均 吸收 系数 (单位 :m '), (c) 


表 5 


Ma =35 热 辐 射 非 耦 合 与 耦合 
位 : MW/m2) 


计算 的 驻 点 热流 值 ( 单 


Table 5 Radiative heat transfers at the stagnation 


point of uncoupled and coupled radiation simulations 


with Ma =35 (unit: MW/m’) 


Heat flux uncoupled coupled thin coupled dom 
qe 2.63 2.51 2.70 
gq; 3,72 2.52 2.92 
4 结论 
本 文 设计 了 辐射 平衡 和 温度 间断 两 个 一 维 热 辐 
射 输 运 模型 算 例 ， 对 比分 析 了 各 类 热 辐射 输 运 求解 


方法 的 性 能 。 同 时 考察 各 方法 在 实际 求解 过 程 中 的 
dt NN tg 
切 平板 近似 和 计算 量 庞大 的 蒙特 卡 
由 离散 坐标 法 、P; 近似 和 光学 厚度 极限 近 0 
同 层次 方法 构建 热 辐 射 输 运 解 算 方法 库 ， 
了 该 方法 库 的 使 
热 辐 射 输 运 ， 能 与 流 场 数值 格式 高 
适合 大 规模 并 行 计算 。 


行 器 高 温 热 


洛 模拟 ， 提 出 了 


本 


同 时 给 出 
原则 。 该 解 算 方法 库 可 准确 求解 
效 耦 合 求解 ， 并 


利用 构建 的 热 辐射 输 运 解 算 方法 库 耦 合 再 入 发 
化 学 非 平衡 流 场 数值 求解 格式 实现 了 


Ma =35 的 球 头 算 例 演示 模拟 。 相应 再 入 飞行 条 件 下 

流 场 与 热 辐射 耦合 效应 显赫 , 热 重 射 对 流 场 具有 ”“ 冷 

却 ” 作 用 ; 流 场 非 耦合 与 耦合 热 鲁 射 计算 的 对 流传 

热 接近 ， 但 非 耦合 计算 的 辐射 热流 远大 于 耦合 计算 

值 。 同 时 ， 耦 合 模拟 预测 的 辐射 传 热 与 对 流传 热 水 

平 相当 。 开 展 相 关 模 拟 研究 时 ， 应 当 考 虑 流 场 的 热 

辐射 耦合 效应 
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